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in 
o 

, The electric and weak electric dipole form factors for heavy fermions are cal- 

culated in the context of the most general two-Higgs-doublet model (2HDM). 
Q H ' We find that the large top mass can produce a significant enhancement of the 

pi ■ electric dipole form factor in the case of the b and c quarks. This effect can 

be used to distinguish between different 2HDM scenarios. 



I. INTRODUCTION 

One of the simplest extensions of the Standard electroweak model (SM) is the so-called 
two-Higgs-doublet model (2HDM), in which the new ingredient is the presence of a second 
doublet of scalar fields. The inclusion of this new fields implies various phenomenological 
consequences, which have led the 2HDM to be subject of analysis during the last two decades. 

We concentrate here on one of the most interesting features concerning the 2HDM, which 
is the presence of different sources of CP violation beyond the standard Sckm phase in the 
quark-mixing matrix. In particular, we analyse the effects of the new parameters of the 
model on the CP-violating electric and weak electric dipole form factors for heavy fermions. 
The interest on these observables has increased in recent years, in view of the ongoing 
activity both in the theoretical and experimental areas SM predictions for CP-odd 

dipole moments are extremely small, and this opens the possibility for one-loop effects 
coming from extended models to show up Specific observables have been proposed 

and studied in the literature and some bounds have already been obtained from 

experimental measurements in e + e~ collisions MM. 
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If done in a completely general way, the addition of a second scalar doublet to the SM 
Lagrangian is problematic: one immediately finds that a general 2HDM model contains 
tree-level flavour-changing neutral currents (FCNC), which are strongly suppressed phe- 
nomenologically. To avoid this problem, it is usual to introduce ad hoc discrete symmetries, 
in such a way that all fermions of a given charge couple to only one of the doublets [Q . This 
can be done in different ways, leading to so-called 2HDM I and II. It is often said that the 
obtained flavour conservation is "natural". The inclusion of discrete symmetries, however, 
is not the only way of preventing the undesired FCNC [pHTTH. In fact, the presence of strong 



hierarchies in the fermion masses and mixing angles seems to be a clear signature of an 
underlying theory of flavour beyond the SM Yukawa couplings. From this point of view, 
it can be also "natural" to expect that the suppression of FCNC observed at low energies 
could be explained in the context of this by now unknown theory. On the other hand, 
whereas the phenomenological constraints on FCNC are very stringent for processes which 
involve the first family of quarks and leptons, this is not the case if one considers only the 
mixing between the second and third fermion families. One possibility is to assume that the 
suppression of FCNC is related to the masses of the involved fermions, as has been proposed 
by several authors in the literature [|K|[iTp • 

Here, instead of choosing a particular Ansatz to enforce the suppression of tree-level 
FCNC, we will consider a completely general 2HDM, using a convenient parametrisation 
to take into account the existing phenomenological constraints. As has been pointed out 
recently [|12|], the various sources of CP violation can be classified into four classes: 



• CP violation in charged and neutral flavour-conserving scalar exchange 

• CP violation in neutral flavour-changing scalar exchange 

• CP violation in the neutral scalar mixing matrix 

• CP violation in charged gauge boson exchange (the usual CKM mechanism) 

It is clear that particular 2HDMs show in general different patterns for these CP violation 
sources. In this paper we analyse and compute the 2HDM predictions for the flavour- 
diagonal CP-odd couplings of heavy fermions, both quarks and leptons, to the neutral 
gauge bosons 7 and Z. For on-shell fermions and gauge bosons, the corresponding ff~f and 
ffZ form factors are known as the electric dipole moment (EDM) and weak electric dipole 
moment (WEDM) of the fermion / respectively. The presence of a nonvanishing dipole 
moment of this kind is a signal of time reversal symmetry violation, and in our framework, 
of CP violation. In general, the form factors are gauge invariant quantities — and can be 
contrasted with experiment — only when the external fermions and gauge bosons are on 
the mass shell. However, it can be seen that in the 2HDM the one-loop predictions for the 
electric and weak electric dipole form factors are still gauge invariant when the gauge bosons 
are off shell. 

We will present the analysis for the CP-violating dipole form factors in the general 2HDM 
case, and then apply the results to models which include discrete symmetries and models 
in which the magnitude of the FCNC is related to the masses of the involved fermions. As 
stated, the fermion electric and weak electric dipole form factors are nonzero only at high 
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orders in the SM, whereas 2HDM give rise in general to nonvanishing contributions at the 
level of one loop. Therefore, they represent good candidates for an observation of CP-odd 
effects arising from an extended scalar sector. In the case of heavy fermions, the effects are 
particularly important, since the new (scalar-mediated) contributions are proportional to 
nonnegative powers of the fermion masses. 

The paper is organized as follows: in section II the notation and the general 2HDM 
considered in this work are presented. The analytical and numerical results for the CP- 
violating form factors are given in section III, while section IV contains our conclusions. 
In the appendix we quote some explicit expressions for the Feynman integrals used in our 
analysis. 



II. MODEL 



As stated, we consider here a completely general 2HDM, allowing in principle for the 



presence of tree-level FCNC. We adopt in the following the notation introduced in Ref. Hl2 
where the most general Higgs potential is parametrised as 



V(01, 02 ) = - ^0202 - f>120102 + h.C. 

+Ai 

1 



f 10l) 2 + A 2 (0 2 2 ) 2 + A 3 (0l01 4fc) + A4(0l02 4<f>l) 



+- 



H4h) 2 + h.c. + (\ 6 4<1>1 + h<t>2<t>2) (0102) + h.c. 



and the neutral scalars acquire vacuum expectation values 



V2 



e i5 cos(3 : 



V2 



sin (3 



(1) 



(2) 



In order to write the scalar-fermion couplings, it is convenient to introduce a new basis for 
the scalars, namely 



^ = e i(5 (cos/3$i + sm/3$ 2 ) 
i 2 = sin/?*! - cos/5$ 2 



with 



+ H° + tG°)/V2 



$5 



(R + i 



(3) 



(4) 



It is easy to see that are physical charged scalar particles, while and G° are the 
Goldstone bosons corresponding to the spontaneous gauge symmetry breakdown. The re- 
maining neutral scalars H°, R and / are not in general mass eigenstates. In terms of these 
fields, the scalar-fermion couplings can be written as 



(V2G F ) l/2 {LM + V2LW) 



(5) 



with 
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L (nt) = ((j lM u Ur + D L M d D R + L L M l L R ) H° + 

(D L T d D R + L L T l L R ) (R + il) + U L T U U R (R - il) + h.c. , 
L {ch) = U L V CKM T d D R H + - D L V ] CKM Y U U R H- + N L V l T l L R H + + h.c. , (6) 

where we have used the definitions U = (u,c,t) T , D = (d,s,b) T , L = (e,/i, r) T and iV = 
(u e , z/ M , v T ) T . The couplings of the Goldstone bosons G° and G ± are the standard ones. As we 
have checked that they do not contribute to the form factors we are interested in, we have not 
written them explicitly in the Lagrangian (|B|). As usual, M* , with / = u,d,l stand for the 
quark and lepton diagonal mass matrices and Vqkm is the Cabibbo-Kobayashi-Maskawa 
matrix (the subindex CKM will be omitted in the following to simplify the notation), 
whereas r*, / = u, d, I are arbitrary 3x3 complex matrices that arise from the extended 
Yukawa couplings. It is useful to distinguish between the diagonal and nondiagonal elements 
of r*, defining j[2| 



(V 



(7) 



where £^ and fi{j are in general complex numbers. The parameters are responsible for 
the tree-level flavour-changing neutral currents. 

Nonzero phases in £f. and represent new sources of CP violation beyond the standard 
Sqkm phase in V. In addition, a further source of CP violation arises from the neutral scalar 
mixing: in the limit where CP is conserved, the CP-even states H° and R do not mix with /, 
which is CP-odd; however, the nonhermitian terms in the Higgs potential can induce either 
explicit or spontaneous CP violation [the latter, arising from the phase 5 in (0)]. Then, 
in general, one expects the neutral scalars to become mixed. The physical neutral mass 
eigenstates H® (i = 1, 2, 3) can be written as 

H?= £ Si S, (8) 

S=H°,R,I 

being O an orthogonal (real) matrix. It is clear that in general H® are not eigenstates of 
CP. 



With the introduction of discrete symmetries to prevent FCNC, the above general struc- 
ture becomes simplified. By requiring the Yukawa couplings to be invariant under the 
changes 

01 -> -01 , 02 -> 02 , 

Dm — > — D Ri , L Ri — > — L Ri , (9) 

together with Uju — > —U R i (Uju — ► U R i), one obtains the so-called 2HDM I (II). The 3x3 
T matrices in (|6|) are then given by 

r*'^M*<, r^[^JJ MU J^j;) . ( 10 ) 
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If the Higgs potential is also invariant under the transformations all CP-violating terms 
in ([!]) turn out to be forbidden. However, one can allow for a soft breakdown of the discrete 
symmetry, retaining the CP violation only through the coupling with //f 2 in ([ft). The electric 
and weak electric dipole form factors for the top quark have been analyzed within this scheme 
in Ref. j|. Notice that, in this case, the parameters in (0) satisfy Im^ = ^ = 0, so that 
the only source of CP violation beyond the SM is the mixing of CP-even and CP-odd fields 
in @. 

As commented above, some models propose specific relations between the magnitude of 
the FCNC and the masses of the involved fermions. One usual Ansatz is that proposed by 



Cheng and Sher [I0|, in which the matrix elements of T u ' are governed by the order of 



magnitude of the fermion masses, obeying 

r"/ = \ i j^m i m j , (11) 

with Xij not far from unity. In general, if the lightest scalar masses are assumed to be in the 
region of a few hundreds GeV, bounds from AF = 2 processes (F = S, C, B) constrain the 
couplings \ s d, X s b and X uc to be < 0.1 Nevertheless, present experimental information 



does not provide such kind of constraints for X ct , which in principle is allowed to be 0(1). 
Phenomenological consequences of having large c-t flavour changing couplings have been 
studied recently by several authors |L4| . We show below that the assumption of the Cheng- 



Sher Ansatz of Eq.flTTD with X ct ~ 1 leads to a significant enhancement in the electric dipole 
form factors of the c quark. 



III. ANALYTICAL AND NUMERICAL RESULTS FOR CP VIOLATING 

DIPOLE FORM FACTORS 

The most general Lorentz invariant matrix element for the interaction of a gauge boson 
B with two on-shell fermions /, / can be written as 

(f(p-)f(pM(0)\B(q))= i e S/ (p_) [if'Vh^ + ^V)7„75 

+ (>#'V) + *&V)tO % + (i#V) + Fl'Vhs) <V<f ] v f (p + ) , (12) 

B f 

where q = p + + p_, e is the proton charge and the coefficients Fj , so-called form factors, 
are in general functions of q 2 . At the tree-level, only the vector and axial vector form factors 
can be different from zero in a gauge theory. 

B f B f 

The last two coefficients, F M and F E , are known as magnetic (weak magnetic) and 
electric (weak electric) dipole form factors when one considers the coupling with the gauge 
boson B — 7 (B — Z). Both F^ and F§^ are chirality-flipping quantities. Here we are 

B f 

interested in particular in F E ' J , which is CP-odd and vanishing small in the SM. We remark 
that within the SM all the above form factors Ff ' (q 2 ) are gauge independent only when the 
gauge boson B is on shell. In this case F E is called the electric (B = 7) or weak electric 
(B = Z) dipole moment of the fermion /. 

We will analyse the values of F E ^(q 2 ) (B = j,Z) for heavy fermions in the 2HDM. 
Within these models one gets in general contributions already at the one-loop level, and 
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the form factors are still gauge invariant when the gauge bosons are off-shell. Since these 
contributions are due to the exchange of neutral and charged Higgs bosons (which carry the 
CP violation effects beyond the SM), and the Higgs-fermion couplings are proportional to 
the corresponding fermion masses, only heavy fermions are expected to provide significant 
effects. We will concentrate in particular on the electric and weak electric form factors for 
the r lepton and the t, b and c quarks. 

In the most general 2HDM, the relevant diagrams that contribute to F E at one loop 
are shown in Fig. 1 [notice that those in Fig. 1 (d) and (e) only contribute to the Zff form 
factor]. Let us begin by quoting the results for F^^(q 2 ) being / an up-like quark. The 
contributions from the diagrams of Fig. 1 (a)-(e) are given by 



(a) 



2^2 d 



sin 9w cos 9w 



9v 



(b) 



(d) 



x E ™/'E mpHOjSffilmCfO^-Re^O^-Reir^^ORjOrj 

f'=u,c,t j=l 
1 



+ -\m(T u ff ,T u flf ) 1(0 



2 — • ./T'./ 7 
2y/2G F 



(O " 2 



I® (mf,mfi, m H o , q 2 
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sin 8w cos 9w 

x £ m f/ Im[(r^V) ff iVT% l )]l( 1 \m f ,m f/ ,m H± ,q 2 ) 



f'=d,s,b 
2^2 G F 



g 



9v m z m f 



sin 9w cos 9w 

x Yl m f/ lm[(T^V) ff (Vr d y ff/ )}l^( mf ,m ft , mH± ,q 2 

f'=d,s,b 

2V2G F 

sin 9w cos 9w 

3 

x J2( lm 0°Rj ~ Re £f°ij) °H°j I m {m f , m H o, q 
3=1 

\f2G F sr- (\ 12 

9dnfl r^O mf E O/Z'l 
Z Sin t7iy COS fyj/ 



(13a) 



(13b) 



(13c) 



(13d) 



71 



f'=u,c,t 



x 



{0 Rj O Rk + OijOi k ) I {W) (m f , m f , m H o , m H o, q 2 ) 



(13e) 



where 



9v 



.1 _ 4 
2 ~ 3 



sin 2 fl 



9 d v 



1 2 

2 + 3 



'-- + -sin 2 fl 



(--+ sin 2 ^) 



and the Feynman integrals are quoted in Appendix A. The gauge invariance of the 
form factors has been explicitly checked. For down-like quarks and leptons, the resulting 
expressions are similar to those in Eqs.(13). In the case of the down quarks these are obtained 
just by replacing 
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T u , g u v <— > T d , g d v 

E — E 

f'=u,c,t f'=d,s,b 

Tj — > -0 7j 

y — ► (14) 

For the charged leptons, the diagrams in Fig. 1 (b) and (c) are zero in the limit of vanishing 
neutrino masses. The remaining contributions can be obtained from those in Eqs. ( |T3D 
through the changes 

9v 

E — E 

f'=u,c,t f'=e,ii,T 

On — > -On (15) 

where g l v = -1/2 + 2 sin 2 9 W . 

Finally, the contributions to the electric dipole form factors are also obtained easily 

Z f 

from the corresponding expressions for F E . In this case the rule simply consists in the 
replacements 

U{ ' ^ll ■ n 1 (16) 



sin 9w cos 9w sin 9w cos 9w 

for f — u,d, I. As stated, the diagrams in Fig. 1 (d) and (e) do not contribute to . 

Now we can make use of these results to evaluate the leading contributions to the electric 
and weak electric dipole form factors for the t, b and c quarks and the r lepton. The final 
expressions can be written as 

3 

f e ' (? ) = E( a i a j + d i ii 

i=l 
3 



W) = E« 



■3 a 3 
i=l 



F|V) = (& + W + E<' 6 7i 

^•V) = (fc + cW+E°J' 6 ^ 

i=l 

i=l 

i f 3 ,T (9 a )=Ea? T aJ 

i=l 
3 

ttiZ.c/ 2\ / / £,c /c ,Z,c c\ 

*V (? ) = E(s- a i %J 

i=l 
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= E («J' C «, C + «r C <) + ( & ' + C ') 7 ' C /?' C (17) 

i=l 

In this parametrisation, the dependence of F E with the unknown quark mass-matrix pa- 
rameters and Higgs-mixing angles has been collected in the factors in Greek letters, whereas 
the coefficients kf , with k = a, a', ... d, contain the global factors including gauge boson 
couplings and fermion and Z masses, plus the Feynman integrals, which depend on q 2 and 
the masses of the Higgs bosons. The letters a, b, c, d identify the diagram from which each 
contribution originates, according to the notation in Fig. 1 and Eqs. (0). We have distin- 
guished with primes the contributions that include tree-level flavour-changing effects. 

In Eqs. ( [T?| ) we have quoted only the dominant terms arising from the expressions (|13|) . 
hence the contributions proportional to light fermion masses have been neglected. In addi- 
tion, the contributions from the diagram in Fig. 1 (e), which are proportional to flavour- 
changing parameters, have been neglected in comparison to flavour-changing terms arising 
from the diagram in Fig. 1 (a) [notice that ( p.3e[ ) vanishes when the matrices are hermi- 
tian]. 

The explicit expressions for the factors in Greek letters can be easily obtained from Eqs. 
(P|). We find 



O H0j [0 Rj hn£f - e Tj ReO] - e Re^O^ + ^lm(^) (0% - 0%) (18a) 

(5 h = -Im(66) (18b) 
7/ = O H0j (0 Rj ImC/ - e Tj ReC/) (18c) 

a' c = - Re(/Wte) R p ld + - Im(/irf/ite) (0 2 Ri - Ol) (18d) 
m c m t 2 m c m t J 3 

P = M \V cb \ 2 hn(U b ) + ^ lm ^ bV - V ^ + (18e) 

where e = +1 for f — u,t and e = — 1 for f = b,r. If we assume that is not very different 
from one for all fermions (as is the case in most models for the quark mass matrices), all 
parameters in Eqs. (|18a ) - (|18c|) are expected to be 0(1). Then, if no accidental cancellations 



occur, the order of magnitude for the electric and weak electric dipole form factors will be 
given by the coefficients kj in Eqs.(|l7]). On the other hand, in the case of the c quark 
we have to deal with contributions proportional to a' c and (3 fC , which contain the flavour- 
changing parameters fi ct , fitc and fi sb . These contributions depend on the Ansatz chosen for 
the quark-mass matrix, and can be very important due to the large top-quark mass. 

To estimate the order of magnitude of the CP-violating form factors in different 2HDM 

B f 

scenarios, we have calculated numerically the values of the coefficients kj for different 
values of q 2 and the Higgs masses. Our results are presented in Tables I-IV. For the b, c, 
and r form factors we have chosen ^fq 1 = 10, 92, 170 and 500 GeV, corresponding to the 
approximate centre-of-mass energies in S-meson factories, LEP1, LEP2 and future e + e~ 
colliders respectively. For the t quark we have taken yfq 2 = 500 and 1800 GeV, the latter 
corresponding to the high-energy pp collider at Fermilab. We have considered neutral and 
charged scalar-boson masses of 100 and 200 GeV. In general, the Feynman integrals are 
expected to be suppressed if the masses of the involved scalars increase, so that the sums 
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in (|TyD should be dominated by the contribution from the lightest Higgs. Notice that in 
the limit where the neutral scalar sector is degenerate in mass, the contributions from the 
diagrams Fig. 1 (a), (d) and (e) to both the electric and weak electric dipole form factors 
vanish owing to the orthogonality of the neutral-scalar mixing matrix O. 

The contributions of the different CP violation sources in the general 2HDM can be easily 
read from Eqs. fljjp . Let us come back to the classification presented in the introduction: the 
CP violation in charged and neutral flavour-conserving scalar exchange is contained in the 
imaginary part of the parameters whereas the CP violation in neutral flavour-changing 
scalar exchange is due to the imaginary part of [iff- The CP violation in the scalar mixing 
matrix arises from the mixing between the CP-odd scalar / and the CP-even fields H° and 
R, that means, the products OnojOij and O^Oij in (|i~8f) . Finally the Sckm CP-violating 
phase in the V matrix only appears together with tree-level flavour-changing parameters in 
Eq. (HD. 

In Table V, we quote the expected orders of magnitude of the electric and weak electric 
form factors at q 2 = m| for the r lepton and the b and c quarks, in both 2HDM I /ll and 
Cheng-Sher-like scenarios. For the 2HDM I and II, which include discrete symmetries to 
prevent FCNC, some of the contributions in (|17D vanish. As stated in the previous section, 
in this case Im£y = \ifp = 0, hence fif = or = /3 = 0, and all the effect arises from 
the CP violation in the Higgs-mixing matrix. On the contrary, in the quark-mixing scheme 
proposed by Cheng and Sher, all terms in Eqs. (|i~7|) have to be considered if the \j parameters 
of Eq. (|ll|) are complex numbers of order one. For the case of the b and c quarks, some of 
the terms that vanish in the 2HDM I and II have relatively large coefficients, proportional to 
the square of the top-quark mass. This shows up in the value of the corresponding electric 
dipole form factors, where the predictions in the Cheng-Sher scheme are about three orders 
of magnitude higher than in the 2HDM I /II, as can be seen in Table V. In the case of the 
weak electric dipole form factors the effect is hidden due to the presence of other important 
contributions proportional to m 2 z . 

For the top quark the CP-violating form factors in the general 2HDM are dominated 
by the contributions of diagrams (a) and (d) in Fig. 1. These arise from flavour-conserving 
neutral Higgs exchange, and do not vanish in general in 2HDM I and II. As stated in the 
previous section, F^ 1 and Fg* have been analyzed previously for these models in Ref. [|J. 
The values for F^' t {q 2 ) arising from Eqs. ( |13a| ) and ( |13d| ) and the equivalent expressions for 
F] ; ,t (q 2 ) are in agreement with the results obtained in that paper. 

The weak electric dipole moments of the r lepton and b quark have been also calculated 
recently within the Minimal Supersymmetric Standard Model (MSSM) [|[. The results are 
similar to those obtained in our general 2HDM scheme: |Re[F^' r (m|)]| < 0.3(12) x KT 21 
ecm, |Re[F^' 6 (m|)]| < 1.4(35) x 10~ 21 ecm. In Ref. the authors also quote the values 
obtained for the top quark form factors at y/q 1 = 500 GeV, which yield approximately 
\F^' t \ ~ l-Fg^l — 10~ 19 ecm. Taking the corresponding coefficients from Table I, this order 
of magnitude agrees with our results. Notice that in our analysis we have not considered 
possible corrections arising from the running of the quark masses (we have used the running 
masses m q in the MS scheme, with /i = m q ). In any case, these effects would not modify 
the quoted orders of magnitude for the values in Tables I-V. 
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IV. CONCLUSIONS 



We compute the CP-violating electric and weak electric dipole form factors for heavy 
fermions in the framework of a completely general 2HDM. Despite of being one of the 
simplest extensions of the SM, this model contains interesting new features, such as the 
presence of various sources of CP violation beyond the standard CKM mechanism. The CP- 
violating dipole form factors are vanishingly small in the SM, thus they are good candidates 
to provide observable signals of new physics. In the 2HDM, at one loop, they are found to 
be finite and gauge invariant quantities, even when the involved 7 or Z bosons are off-shell. 

The effect of the different sources of CP violation on the form factors for the c, b and 
t quarks and the r lepton is shown in Eqs. ((T~8|) . In particular, it is seen that some of the 
contributions vanish for the so-called 2HDM I and II, which include discrete symmetries to 
eliminate undesired FCNC. In these models the only remaining terms are those involving the 
mixing between CP-even and CP-odd neutral scalars. In the general case, however, these 
symmetries may be not present. If no accidental cancellations occur, this would imply an 
enhancement of three orders of magnitude in the electric dipole form factor of the b quark 
with respect to the prediction of 2HDM I and II for energies in the GeV to TeV range. On 
the other hand, we find that in the case of the c quark the electric dipole form factor is 
strongly dependent on the presence of c-t flavour-changing effects. Assuming an up-quark 
mass matrix of the type proposed by Cheng and Sher, with \ uc ~ 0(1), the values for F]/ c 
for energies from 10 to 500 GeV can be two to four orders of magnitude larger than those 
obtained in 2HDM I or II. We conclude from this analysis that the study of CP-violating 
dipole form factors can be useful to get information on the flavour mixing, and offers an 
interesting possibility to distinguish between these and other possible 2HDM scenarios. 
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APPENDIX A: FEYNMAN INTEGRALS 



The integrals /W introduced in Eq. flT3| ) are defined as follows: 




(Al) 



y) = ml(l-x-y) + {m 2 q/ - m 2 q ) (x + y) + m\ (x + yf 



— sxy 



Im[/ (I) (m g ,mg/,m0,s)] 



1 /V 




1 ( P 2 g ~P/ 




16ns [3 2 q 



PA' V 4 
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x log 
1 

16tt 2 



' s (P q + P q ,f + Ami 
s{(3 q -(3 ql y + Aml 
i 

o 



Q(s- Ami) 



P.V 



dx 



o 



!-* 1 - 2x 
rf ^T7 V 



f 2 {x, y) = m 2 ql (l-x-y) + (mj - m 2 .) (x + y) + m 2 , (x + y) 2 - s x y 



Im[/ (II) (m q ,m q >,m^, s)} 



Ins (3? 



J_(Pl + Pl + ^ 



Re[/( m) (m g ,m^ S )] 



x log J 
1 

_ 16tt 2 



,s(/3 g -/^) 2 + 4m 2 ,, 



RV. / rfx 



6( S -4m 2 ) 
2/ 



f 3 (x,y) = m 2 q (1 - x - y) 2 + m 2 ^x + m 2 z y - s x y 



Im[/ (III) (m„m^s) 



1 



16ns \/3 2 2(3 q 
c - b z /3, 



' c-P , 



Re[/ (IV) (m 9 , m q/ , m<y, s)] 



x log 
1 

~16tt 2 



c + b z (3 q 



P.V. 



2 2 2 

s 

\2l 



c/.r 



0[s - (m z + m^) } 

(x-y)(l-x-y) 



l-x 



dy- 



fa(x,y) 



h(x,y) = {ml - m 2 ) (1 - x - y) + m 2 (1 - x - yf + mix + ml,y - sxy 



Im[/ (IV) (mq, nty, m , 77y , s)] 



1 (m^-mp) 



16tt 



2V, d i 



d - bp P Q 



d+ bs'Pa 



Q[s - (m<p + 



with 



/? a = (l-4m 2 / S ) 1 / 2 , a = q,q',<f> 

b a = [(1 - (m + mfl ) 2 /sf/ 2 [1 - (m* - ffl„) 2 /s] 1/2 , a = Z, 0' 
c = 1- (mj + m|)/s 
d = 1 - (mj + mj,)/s - 2 (m 2 - m 2 ,)/s 



(A2) 
(A3) 



(A4) 
(A5) 



(A6) 
(A7) 



(A8) 



(A9) 



In these expressions we assume s > Am 2 , which is valid for all the numerical estimations 
that are presented in this paper. The analytical results in (A~2) and (|A6|) agree with those 
presented in Ref. || in the limit (3 q i = Pq- 
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TABLES 



"Z t "7 t 

TABLE I. Dominant coefficients for F E ' (s) and (s) for different values of s and the Higgs 
mass ran- All values are in units of ecm. 





= 500 GeV 
= 100 GeV m H = 200 GeV 


y/s = 1800 GeV 
m H = 100 GeV m H = 200 GeV 


a z,t 


(-0.1 + 1.4*) x 10" 19 
(1.0 - 2.8*) x 10" 20 


(0.2 + l.Oi) x 10~ 19 
(0.8 - 3.5*) x 10" 20 


(-1.2 + l.Oi) x 10" 20 
(3.1 - 2.6i) x 10" 21 


(-1.0 + l.Oi) x 10" 20 
(3.0 - 2.6i) x 10" 21 




(-0.3 + 4.1*) x 10" 19 


(0.6 + 2.9*) x 10~ 19 


(-3.4 + 2.9i) x 10" 20 


(-2.9 + 2.7*) x 10" 20 



TABLE II. Dominant coefficients for F§ fi (s) and F% b (s) for m# = 100 and 200 GeV and 
different values of s. All values are in units of ecm. 





y/i = 10 GeV 


\/s = mz 


s/s = 170 GeV 


y/a = 500 GeV 




m H = 100 GeV 


(b + c) z > b 

d z > b 


5.0 x 10~ 21 
3.6 x 10~ 21 


5.1 x 10~ 21 
4.0 x 10~ 21 


5.7 x 10~ 21 
6.3 x 10~ 21 


(2.8 + 5.9i) x 10" 21 
(-1.1 + 1.8i) x 10~ 21 


a'' b 

(b + c)^ b 


(-0.8-0.2i) x 10~ 22 
1.0 x 10~ 20 


(-0.7-2.6i) x 10~ 23 
1.0 x 10~ 20 


(0.2 - 1.4i) x 10" 23 
1.1 x 10- 20 


(2.1 - 2.7*) x 10" 24 
(0.5 + 1.3*) x 10" 20 




m H = 200 GeV 


(b + c) z > b 
d z > b 


3.0 x 10~ 21 
2.0 x 10~ 21 


3.1 x 10~ 21 
2.1 x 10~ 21 


3.3 x 10~ 21 

2.4 x 10~ 21 


(3.4 + 3.2*) x 10" 21 
(-0.5 + 2.5*) x 10~ 21 


a^' b 

(b + cp' b 


(-2.3 -0.4i) x 10~ 23 
6.5 x 10~ 21 


(-0.6 - 0.9i) x 10~ 23 
6.6 x 10~ 21 


(-0.2 - 0.7i) x 10~ 23 
7.0 x 10~ 21 


(0.8 - 2.2i) x 10" 24 
(6.5 + 7.8*) x 10~ 21 



TABLE III. Dominant coefficients for F^ T (s) and F£ T (s) for m H = 100 and 200 GeV and 



different values of s. All values are in units of ecm. 





s/s = 10 GeV 


\fs = m z 


s/s = 170 GeV 


y/Ji = 500 GeV 




m H = 100 GeV 


d z > T 


1.6 x 10~ 22 


1.8 x 10~ 22 


2.9 x 10~ 22 


(-0.5 + 0.9i) x 10~ 22 




(-1.2 - 0.7i) x 10~ 23 


(-l.0-4.8i) x 10~ 24 


(0.4 - 2.7*) x 10" 24 


(0.4 - 0.5*) x 10~ 24 




m H = 200 GeV 


d z > T 


0.9 x 10~ 22 


0.9 x 10~ 22 


1.1 x 10~ 22 


(—0.2 + 1.1*) x 10~ 22 


a"t' T 


(-0.4-0.2i) x 10" 23 


(-1.1 - 1.6*) x 10" 24 


(0.4 - 1.2*) x 10" 24 


(0.1 - 0.4i) x 10" 25 
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TABLE IV. Dominant coefficients for F^ c (s) and F% c (s) for m H = 100 and 200 GeV and 



different values of s. All values are in units of ecm. 





V« = 10 GeV 


V 7 * = ™Z 


V« = 170 GeV 


yfe = 500 GeV 




m# = 100 GeV 


/Z,c 

a 

AZ.c 

a ' 


0.6 x 10" 21 
-0.7 x 10~ 21 


0.6 x 10~ 21 
-0.7 x 10~ 21 


0.7 x 10~ 21 
-1.2 x 10~ 21 


(0.2 + 1.0s) x 10" 21 
(2.0 - 3.5s) x 10" 22 


a'' c 
a n ' c 

{V + c') 7,c 


(5.9 + 2.8i) x 10~ 24 

1.7 x 10~ 21 
(5.3 + l.Oi) x 10" 24 


(0.5 + 1.9*) x 10- 24 

1.8 x 10~ 21 
(1.0 + 1.6?) x 10" 24 


(-0.1 + 1.1?) x 10~ 24 

1.9 x 10~ 21 
(0.6 + 0.9?) x 10~ 24 


(-1.6 + 2.0?) x 10~ 25 
(.6 + 2.9?) x 10~ 21 
(0.5 + 7.5?) x 10" 25 




m H = 200 GeV 


a /Z ' c 
d Z,c 


4.4 x 10~ 22 
-3.7 x 10~ 22 


4.6 x 10~ 22 
-3.8 x 10" 22 


4.9 x 10~ 22 
-4.4 x 10~ 22 


(2.5 + 7.2?) x 10" 22 
(1.0 - 4.6?) x 10" 22 


a<' c 
a' 7 ' c 

{b' + c') 7,c 


(1.6 + 0.7?) x 10" 24 

1.3 x 10~ 21 
(1.6 + 0.3i) x 10" 24 


(4.3 + 6.4?) x 10~ 25 

1.3 x 10~ 21 
(0.5 + 0.5?) x 10~ 24 


(1.4 + 5.0?) x 10- 25 

1.4 x 10~ 21 
(2.7 + 4.1?) x 10~ 25 


(-0.6 + 1.6*) x 10~ 25 
(0.7 + 2.1?) x 10- 21 
(-0.2 + 1.4?) x 10" 25 



TABLE V. Expected order of magnitude for electric and weak electric form factors at q 2 = m| 
for r lepton and b and c quarks in 2HDM I /ll and Cheng-Sher-like scenarios. Values are in units 
of ecm. 









F Z E '\m\) 


F'Aml) 


F Z E ' c {m\) 


Fr^i) 


2HDM I/II 


1Q -22 


1Q -24 


lO" 21 


1Q -23 


lO" 21 


1Q -24 


Cheng-Sher 


1Q -22 


1Q -24 


lO^-lO -21 


1Q -20 


10 -21. 10 -21| Act |2 


10" 21 |A ct | 2 
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